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To visualize a higher dimensional object it is convenient to consider its two-dimensional 
cross-sections. The set of quantum states for a three level system has eight dimen- 
sions. We supplement a recent paper by Goyal et al by considering the set of all 
possible two-dimensional cross-sections of the qutrit. Each such cross-section is 
bounded by a plane cubic curve. 
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1 Introduction 



With the rise of quantum information theory, a minor issue but intriguing issue has 
attracted some attention: What does the set of quantum states actually look like? 

For a qubit the set of quantum states is a ball, but the case of a qutrit is 
already complicated [T|. It is an eight dimensional convex body. A standard way to 
visualize higher dimensional sets is to consider two-dimensional cross-sections, and 
preferably an intelligently organized Grand Tour of such 2-sections [2j. Of course 
this idea has been pursued for the qutrit, typically by considering the finite set 
of cross-sections generated by the Gell-Mann matrices [3] - j8|. An elegant and 
comprehensive treatment is due to Goyal et al [§]. However, by concentrating on 
cross-sections generated by the Gell-Mann matrices one obtains only 28 2-sections 
altogether, or 56 three-dimensional 3-sections. This is quite far from a Grand Tour. 
On the other hand we are looking at a fairly symmetric eight-dimensional body, left 
invariant by an SU(3) subgroup of the rotation group 50(8), and this cuts down 
the size of the problem. Our purpose here is to supplement the work by Goyal et 
al. by making this idea precise. In section [4] this will enable us to see how their 
selection samples the full set of 2-sections, because we will be able to visualize the 
set of all unitarily inequivalent two-dimensional cross-sections of the state space of 
the qutrit. In section [5] we will consider the shape of a general 2-section, and the 
plane cubic curves that bound them. 

Because of the duality between sections and projections of a self-dual set, a 
description of two-dimensional cross-sections can be translated to results on the 
shapes of the numerical ranges of complex matrices [10] . The numerical range is an 
interesting tool which has found some applications in quantum information theory 
[11]. This gives some special interest to the study of 2-sections, as opposed to 
3-sections. We discuss this topic in section [6] 



2 Preliminaries 

A 2-section of an object in an eight dimensional vector space is its intersection with 
a 2-plane through the origin — which will be placed at the maximally mixed qutrit 
state in our case. The modeling linear subspace for such a plane is a 2-dimensional 
subspace V of the real vector space of traceless 3x3 hermitian matrices equipped 
with the Hilbert-Schmidt scalar product 

M 1 -M 2 = -TrMM . (1) 

The set of quantum states is the convex set of traceless hermitian matrices M such 
that the density matrix p — 1/3 + Mis positive. It is inscribed in a minimal sphere 
of radius R ou t, and contains a sphere of maximal radius R[ n , where 

-Rout = ^/^' = 2"^ out ' ^ 

We need to know the dimension of the set of all 2-sections, or equivalently of 
the set of all 2-dimensional subspaces of a real 8-dimensional space. This set is also 
known as the Grassmannian Gr(2, 6). Its dimension is easily found: 

Fact 2.1. The set of orthonormal bases in M. n is isomorphic to the orthogonal group 
Q(n). 
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Fact 2.2. The dimension of the group 0(n) is — '-. 

Fact 2.3. The set of 2-dimensional subspaces of an 8-dimensional real vector spaces 
has 12 dimensions. 

Proof. : Any 2-dimensional subspace can be obtained by choosing an orthonormal 
basis such that its first two vectors span the subspace, and the remaining six its 
orthogonal complement. But we obtain the same subspace if we change the bases 
within the subspace and within its complement. Hence the dimension of the set of 
2-planes equals the dimension of 0(8) minus the dimension of the subgroup 0(2) x 
0(6). So the answer is 8 • 7/2 - 2 • 1/2 - 6 • 5/2 = 12. □ 
A general density matrix is given by 

p = i/ 3 + M , (3) 

where .Z3/3 is the maximally mixed state and M is a traceless hermitian matrix 
chosen such that all eigenvalues of p are non-negative. The group SU(3) acts on 
our vector space through M -»■ UMW. This is an 8-dimensional subgroup of the 
rotation group 50(8). Since the action of SU(3) leaves the set of density matri- 
ces invariant all sections with subspaces related by this action will be regarded as 
equivalent. Hence we have arrived at: 

Fact 2.4. The dimension of the set of inequivalent 2-sections of the set of states of 
the qutrit has only 12 — 8 = 4 dimensions. 

Equivalent 2-sections have the same shape, but the converse does not hold [9j. 

3 Representatives of 2-sections 

Let V be a 2-dimensional subspace of the real vector space of traceless hermitian 
matrices. We begin with a simple observation: 

Lemma 1. The subspace V contains no elements of rank 1 and at least one element 
of rank 2 . 

Proof. Assume that the subspace V is spanned by elements A and B. A non-zero 
element of V is either of rank 2 or of rank 3 because is traceless. Then assume that 
A is of rank 3 (otherwise the proof is complete). Consider the determinant of a 
linear combination: 

det(XA + B) = 

Because A and B are hermitian, this is a real polynomial of the order 3. It has at 
least one real root. If A is this root the combination is of rank two. □ 

Since we are interested in picturing the set of 2-sections up to unitary equivalence, 
we use Lemma 1 to introduce a standard form for a basis in the 2-plane. We choose 
the first basis vector A to have rank 2, and then we choose the basis of C 3 so that A 
is diagonal. This fixes the basis in C 3 up to diagonal unitary transformations, and 
we use the remaining freedom to adjust the phases of the orthogonal basis vector 
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B. In this way we arrive at the following basis for the 2-plane: 
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where we may assume that a, b,c,k> 0. We also define d = \/3k. The matrices A 
and B form an orthonormal basis if and only if 

3k 2 + a 2 + b 2 + c 2 = a 2 + b 2 + c 2 + d 2 = 1 . (6) 

The set of possible tuples (a, b, c, d) form 1/16 of a three dimensional sphere, and so 
it is homeomorphic to a three dimensional simplex. The set of matrices B contains 
one representative of each unitary equivalence class of 2-sections of the set of qutrit 
states. We can think of it as a Cartesian product of a simplex and a circle, except 
that the circles shrink to points at three of the four faces of the simplex. 

Unfortunately some ambiguities remain because there may be more than one 
matrix of rank 2 in the 2-plane, and one can then perform a unitary transformation 
so that another matrix takes our standard form A. In particular given a matrix of 
rank 2 its negative also has rank 2. If we leave the origin in the direction of a rank 2 
matrix we will hit the boundary of the set of states at a density matrix of spectrum 
(2/3, 1/3, 0). Such a density matrix lies on a sphere of radius R 2 , where 

Rin < i?2 = < -Rout ■ (7) 

Goyal et al [S] refer to this as the self-dual sphere, because a 2-section containing 
a point on this sphere will also contain its antipodal point. But this contradicts 
another and more common use of the word self-dual, to be discussed in section 
|6j Anyway the point is that the boundary of the set of quantum states intersects 
the sphere of radius R 2 in antipodal points. Now a permutation of the first and 
the second vector of the basis of the Hilbert space will change the sign of A and 
exchange the parameters b, c in B. Thus we are led to identify pairs of points related 
to each other by reflection with respect to the surface b = c. We reduce then the 
set of representatives to the halfsimplex HS of matrices B where b > c (times the 
circle, which shrinks to a point at two of the faces of the halfsimplex). 

To see if there are further ambiguities we solve the equation det(XA + B) = 0. 
By definition of A one solution is for A = oo. The other solutions are the roots of 

2k ■ X 2 + \{b 2 - c 2 ) + 2abc cos - 2A; 3 + k{2a 2 - b 2 - c 2 ) . (8) 

First consider the situation when all elements of the subspace have rank 2. It 
happens if and only if k — 0, b = c and abc cos <fi = 0. The resulting 2-section is 
always a circular disk of radius R 2 . As noted by Goyal et al (in a special case) these 
2-sections are unitarily inequivalent even though they have the same shape. To see 
this, observe that 

\ti{AB - BA) 2 = Aa 2 + b 2 + c 2 = 3(a 2 - k 2 ) + 1 . (9) 
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Once we have restricted ourselves to k = it follows that the value of a cannot be 
changed by changing the basis in the 2-plane V, so no ambiguities arise in this case. 

If the discriminant of ^ is greater than or equal zero we have a discrete ambi- 
guity in the choice of our matrix A, and this does give rise to discrete ambiguities 
in our parametrization. In particular one can check that 



(k,a,b,c) = K),O,cos^,sin0 and (k, a, b, c, 0) = (^—^—, ~^7|~> 

(10) 

correspond to equivalent 2-sections. In general the discriminant of eq. ^ vanishes 
if either k = or a complicated ^-dependent condition holds. We decided to ignore 
this difficulty 

We get the following: 

Theorem 1. Any 2- dimensional section of the space of traceless hermitian matrices 
in C 3 can be represented after appropriate change of basis o/C 3 as a subspace spanned 
by matrices of the form given in eqs. Q^j, where a 2 + b 2 + c 2 + 3k 2 — 1, b > c 
and cf) is arbitrary. Topologically this set forms a simplex times a circle, with the 
circles shrinking to points at two of the faces of the simplex (when abc = 0). The 
parametrization determines the set of unitarily equivalent 2-sections uniquely except 
for discrete ambiguities that occur if the 2-section contains exactly 4 or exactly 6 
traceless matrices of rank 2. 




Figure 1: Halfsimplex. All points not con- 
tained in the dashed area (with its bound- 
ary) has an additional degree of freedom. 





Figure 2: Halfsimplex deformed and 
represented in two dimensions. All 
points from the dashed area are marked 
as red. 



The halfsimplex looks as in figure [T] All points not contained in the dashed 
area (with its boundary) have an additional degree of freedom — the angle 0. To 
visualize this 4-dimensional set, first deform the half-simplex homeomorphically (not 
diffeomorphically) such that all dashed area lie on one 2-dimensional surface, and 
the surface b = c is orthogonal to it. Now reduce the dimension of such a deformed 
simplex as in figure [2} Next we rotate in four dimensions using rotations that 
leave the surface a = V c = invariant. In the reduced picture this is related 
to rotation around the axis y. After rotation we get a half of a closed ball. The 
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topology of resulting set is not interesting, however keep in mind, that to get one- 
to-one correspondence with the set of sections, some pairs of points should be glued 
together. 



4 Sections spanned by Gell-Mann matrices 



Goyal et al [9j consider 2-sections spanned by pairs of the eight Gell-Mann matrices, 
and illustrate them beautifully. The Gell-Mann matrices form an orthonormal basis 
of the vector space of traceless hermitian matrices, and are defined by 



i=l 



xs/y/S + xs x\ — ix2 x^ — ix^ 
Xi + ix 2 x 8 /\/3 — Xs x 6 — ix 7 . (11) 
Xi + ixs xq + ixi —2xs/\/3_ 

Matrices Ai, . . . , A7 are of rank two and are unitarily equivalent to our matrix A = A3. 

To check how a section spanned by a pair of Gell-Mann matrices is represented 
we bring one of them to the standard form A3 and call it A. The same operation 
brings then the other matrix to the form B, which is the representative of the section. 
Calculating it one gets (and we refer to Goyal et al. jS] for illustrations): 

• For pairs 12, 13, 23, 45, 67 the representative is 

1 

1 


and the shape of intersection is circular disk. As a side remark we observe 
that the 4-section given by the quartet 1245 (say) is a round ball. 

For pairs 14, 15, 16, 17, 24, 25, 26, 27, 46, 47, 56, 57 the representative is 



1 

7i 







1 1 



and the shape of the intersection is a circular disk, unitarily inequivalent to 
the above. 

For pairs 34, 35, 36, 37 the representative is 



and the shape of the intersection is a parabola. Note that the discrete ambi- 
guity in our parametrization turns up here. 
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For pairs 18, 28, 38 the representative is 
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7! 



1 
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and the shape of the intersection is an equilateral triangle. 
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• For pairs 48, 58, 68, 78 the representative is 

13 " 
3 1 
0-2 

and the shape of the intersection is an ellipse. 

All these cases have cos (p — and sit at the edges of the halfsimplex, as shown 
(together with their shapes) in figure [3] 




Figure 3: The shapes of special sections and their positions in the halfsimplex tp = 0. 
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5 The shape of a general 2-section 

So far we have sampled only some very special points in the set of all 2-sections. As 
we have seen, up to unitary transformations it is enough to look at density matrices 
of the form p = 1^/3 + xA + yB. The boundary of a 2-section is described by 

3detp = ]--x 2 -y 2 + 3(2abccos<p-k(l-k 2 -3a 2 ))-y 3 + 6k-x 2 y + 3(b 2 -c 2 )-xy 2 = . 
9 

(12) 

The problem of classifying all 2-sections is thereby reduced to the problem of clas- 
sifying plane cubic curves of a somewhat special form. Here we will be concerned 
with two questions: When does the cubic curve factorize into three linear factors 
or into one linear and one quadratic factor? If it does not factorize, when is the 
boundary of the 2-section not smooth? 

The cubic factorizes in the following cases: 

1. b = c A (6 = V acosip = 3k). There are two families of sections. One lies 
on the boundary b = c = of the simplex and the second forms for any G 
(— f, f) a line starting from the point b = c = l/v2 and passing the surface 
b = c. For (p E (|, 3|) the second family is only one point b = c = 1/a/2. The 
2-section is described by 

Q + 2ky - 3x 2 + 3(Ak 2 - l)y 2 ^j > (13) 
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If k — | the boundary is a parabola intersecting a line on the outsphere. 
Otherwise the condition describes a cut hyperbola or a (possibly cut) ellipse, 
depending on the sign of 4k 2 — 1: 

I - 2 "y) ( 3 <« 2 - vt*+ w^t/ - m^T) - 3 * 2 ) - (14) 

If k = this is a circle, and if k — > 4g the cut hyperbola goes into a triangle. 
A section from the second family is always related to an uncut ellipse. 

a = A (b = V c = 0) (if we want to stay in the proper halfsimplex, the only 
possibility is c = 0). These points form an edge of halfsimplex. The cubic 
condition reduces then to: 

i ± x + ky^j Q + 3(k 2 - l)y 2 =f x ± 6A;xy - Jfej/^ > (15) 
If A; = the condition describes a parabola intersected by a line: 

i±z + *2/J Q-VT^ >0 (16) 
If k ^ the condition describes a hyperbola intersected by a line: 

£ ± - + A*) ((^ =F - [ik 2 - l)y ±kx- k/6) 2 + (1 - ~ ^ 

If k — ^ the cut hyperbola degenerates to a triangle. 

Fact 5.1. Assume that the boundary of the section contains a line segment. Then 
the rest of boundary is a conic which connects to the segment in two exposed points. 

We now turn to the second question: if the boundary does not factorize, can it 
fail to be smooth? We will see that this happens if and only if the section contains 
a pure state. 

The cubic curve ( 12 ) will fail to be smooth if and only if 

det p = d x det p = d y det p = . (18) 
A straightforward calculation verifies that this happens if and only if 

x 2 + y 2 = 1 -. (19) 

This is precisely where the boundary of the 2-section touches the outsphere, of radius 
4g, which means that it happens if and only if the 2-section passes through a pure 
state. So we want to find all 2-sections containing a pure state. 
Assume that the matrix 



> 
(17) 



p = J 3 /3 + xA + yB 



- + x + ky yae 1 ^ ybe 1 ^ 



3 

yae~ 1 ^ | — x + ky yce 1 ^ 



ybe' 14, yce- i4 > \ - 2ky 



(20) 



is of rank one. We will use X, Y, Z to denote the diagonal entries. First observe, 
that y ^ 0. Now consider the minor formed by the first and the second column and 
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Figure 4: Points in the simplex ip = containing a pure state 



the first and the third row. We have that Xcye~ 1 ^ = aby 2 . If ^ {0,7r}, then one 
has a = V b = 0. In this case one can remove the phase, and it is enough to 
consider real matrices with e 1 ^ = ±1. 

We have three equations arising from non-main minors: 

Xc = ±aby Yb = ±cay Za = ±bcy 

If one of the numbers a, b, c is zero, then at least one other has to be zero. In this 



case the matrix (20) has a block structure, the determinant factorises, and the pure 
state lies in the points where both factors simultaneously vanish. We have already 
discussed this case. 

Observe that in the first family of sections with factorised boundary both factors 
can vanish simultaneously only if k < 1/2 \/3. Otherwise the boundary is an ellipse 
without a pure state. In the limit case k = 1/2V3, one has an ellipse with one pure 
state on it. 

Consider now the case when abc ^ 0. Consider the equations for the main 
minors: 

XY = a 2 y 2 YZ = c 2 y 2 ZX = b 2 y 2 

One can easily calculate: 

be ab ca 

X = ±-y Y = ±—y Z = ±—y 

a c b 

Using the normalisation of the trace one has y = ( afe )2 + (^ c 2+ ^ 2 . Applying it to 
the above equations one gets: 6kabc = (ab) 2 + (ac) 2 — 2(6c) 2 . Points satisfying this 
equation lie on the surface presented on the figure |4j 

An example of such a section isk = 0,a = b = c= The point representing 
this section lies in the middle of the upper wall of the simplex. 



6 Sections, projections, and numerical ranges 

An important property of the set of quantum states is its self-duality. This gives 
quantum logic its special flavour [T2"j [TBJ HI]- In the vector space of traceless hermi- 
tian n x n matrices the dual (or polar) X* of a set X is defined by 
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X* = {M : 1/n + Tr (MM') > V M' E X} . (21) 

If we recall that any density matrix can be written as p = I n /n + M, where M is 
chosen so that p is positive, it is easily seen that the self-duality of the set of density 
matrices follows from the fact that Trpp' > for any pair of positive matrices p and 
p' . Subsets of the set of density matrices are not self-dual. If we take for example 
the set of separable states, the dual set will be a set of normalized entanglement 
witnesses. The properties of the duality operation are neatly summarized by 



X** = X (22) 

X C Y X* D Y* (23) 

0* = W 1 "- 1 (24) 

(X U Y)* = X* n Y* (25) 

(X n Y)* = conv(X* U Y*) (26) 



where the convex hull appears in the last line. By definition the convex hull of a 
set is the smallest convex set that includes the given set. For a nonsingular linear 
transformation A one has 

(AX)* = (A~ l ) T X*. 

In particular if A is orthogonal the set and its dual transform in the same way. 

The reason why we bring this up here is that a cross-section of a self-dual body 
is dual to a projection onto the linear subspace defining the cross-section j 15 j . Here 
we are interested in cross-sections using planes that pass through the origin (the 
maximally mixed state), and the projection is orthogonal (orthographic). Taking a 
cross-section by means of a plane not passing through the origin will give the dual of 
a perspective projection from a finite point. In mathematics a cross-section of a cone 
of semi-positive definite matrices is called a spectrahedron, and the question what 
kind of convex bodies that can be obtained as projections of spectrahedra arises 
naturally. Once we understand the set of 2-sections of the qutrit we can answer this 
question in our special case [8J. 

There is a further connection to the notion of numerical range of a matrix [TBI 
ITT] . For a quadratic complex matrix A we define a subset of complex plane called 
numerical range of A by 

W(A) = {Tr(Ap) :pen} (27) 

where Q denotes the set of density matrices. Let A = al + B + iC, where B, C are 
hermitian and traceless. Assume first that B, C are orthogonal and of unit norm. 
Then one has 

W(A) = {p B + ip c : peQ} + a 

where ps and pc denote the components of p respectively in the direction of B and 
of C. We get a translated projection of Q onto an afline subspace passing through 
the maximally mixed state. If we abandon the assumption that B, C are normalized 
and orthogonal our set will be a linearly deformed and translated projection of Q. 
The numerical range of 3 x 3 matrices is well understood [17J, and the idea was used 
recently to explore the shadows cast on 2-planes by the set of quantum states |10| . 
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7 Conclusions 



Our intention with this note was to give a parametrization of the set of all 2-sections 
of the qutrit. Theorem 1 is a quite satisfactory answer to this problem. In section 
[4] this enabled us to see at a glance how earlier works — in particular the interesting 
work by Goyal et al [9 J — sample this set. In section [5] we explored how the shape of 
the 2-section changes as we move through the set of all 2-sections, and in section [6] 
we gave a glimpse of a more general context to which the study of cross-sections of 
sets of positive matrices belongs. 
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